Dynamic Nuclear Polarization in Double Quantum Dots 
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We theoretically investigate the controlled dynamic polarization of lattice nuclear spins in GaAs 
double quantum dots containing two electrons. Three regimes of long-term dynamics are identified, 
including the build up of a large difference in the Overhauser fields across the dots, the saturation 
of the nuclear polarization process associated with formation of so-called "dark states," and the 
elimination of the difference field. We show that in the case of unequal dots, build up of difference 
fields generally accompanies the nuclear polarization process, whereas for nearly identical dots, build 
up of difference fields competes with polarization saturation in dark states. The elimination of the 
difference field does not, in general, correspond to a stable steady state of the polarization process. 

PACS numbers: 73.21. La, 76.60.-k, 76.70.Fz, 03.65.Yz 



Understanding the non-equilibrium quantum dynam- 
ics of localized electronic spins interacting with a large 
number of nuclear spins is an important goal in meso- 
scopic physics [JQ- These interactions play a central 
role in spin-based implementations of quantum informa- 
tion science, in that they determine the coherence prop- 
erties of electronic spin quantum bits Q. One of the 
promising systems for realization of spin-based qubits in- 
volves electrically-gated pairs of quantum dots in GaAs, 
with one electron in each quantum dot (Fig. [Tp) 
Hyperfine interactions with lattice nuclear spins are the 
leading mechanism for decoherence of the electron spins, 
and efforts are currently being directed towards under- 
standing these interactions [10Hl5j . with the ultimate 
goal of turning the nuclear spins into a resource by con- 
trolling these interactions [16H19| . Recent experiments 
have successfully demonstrated a wide variety of electron- 
controlled nuclear spin polarization dynamics [l9l- 22|. 
but to date there is no unifying theoretical framework 
in which to understand the experimental results. 



In this Letter we investigate theoretically the process of 
dynamic nuclear polarization (DNP) in two-electron dou- 
ble quantum dots. This process involves the preparation 
of the electronic spins in a singlet state and subsequent 
level crossing between the electronic singlet and triplet 
states with different projection of electronic angular mo- 
mentum (Fig. [T^i) [2fj| | . It is accompanied by nuclear spin 
flips, which polarize the spins of the nuclei inside the 
two dots, producing an effective magnetic (Overhauser) 
field for the electronic spins. Experiments demonstrate 
that DNP strongly modifies the difference between the 
Overhauser fields on the two dots, which is of central im 
portance for control over singlet-triplet qubits ljl, 21 



Detailed understanding of DNP in these systems is both 
of fundamental interest and great practical importance 
for GaAs based electron spin qubits 23-261]. 



In what follows we develop a theoretical framework to 
study the non-equilibrium polarization dynamics of the 
nuclear spin environment. Our approach takes advantage 
of the large effective temperature of the nuclear spins 
and the short time-scale for electron spin evolution to 
coarse grain the electronic system's dynamics, yielding 
a master equation for the nuclear spin degrees of free- 
dom, which we solve in a semiclassical limit. Our key 
results may be understood by first considering three pos- 
sible regimes that result from the DNP process. These 
include i) build-up of an effective difference field, ii) sat- 
uration in so-called "dark states," and iii) preparation of 
nuclear spins in each quantum dot in states that produce 
identical Overhauser fields. 

For example, i) in the case of two dots with unequal 
sizes the growth of an Overhauser difference field D z can 
be understood in the following heuristic picture, which 
is borne out by our analytic and numerical calculations. 
Consider a system with a homogeneous wavefunction in 
the presence of both strong DNP pumping and nuclear 
dephasing. The size difference results in different effec- 
tive hyperfine interactions gt( r ) 011 the left(right) dot. We 
find that the nuclear spins have nearly equal spin flip 
rates on the two dots, so that the build up of the total 
Overhauser field S z is proportional to ge + g r , while the 
build up of D z is proportional to gg — g r . Thus, D z tends 
to grow with S z such that D z /S z —> (ge, — g r )/{gt + gr)- 
On the other hand, ii) when the dots are identical, or 
nearly so, we find a second regime at strong pumping, 
where D z does not grow and the polarization process 
shuts down the growth o f S? . by driving the difference 
field towards a dark state [27], with D x = D y = 0. Such 
states are of interest for use as long-lived quantum mem- 
ory. Finally, iii) electronic and nuclear degrees of freedom 
can be completely decoupled if two electrons are initially 
prepared in the singlet state, while the nuclear spins are 
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FIG. 1: a) Two-electron energy levels as a function of detun- 
ing e between (1, 1) and (0, 2) singlet states. The DNP cycle 
is illustrated by arrows, b) A double quantum dot with two 
electrons interacting with a large number of lattice nuclear 
spins, c) Electronic energy level diagram with transitions 
from s to triplet states T+, To, T- driven by Overhauser fields 
D- , D z , D+ respectively (gray arrows) and energies from ex- 
ternal field Bext and exchange splitting J between s and To 
(black arrows). When D± — electron- nuclear flip-flops are 
prevented, and when D — 0, electrons and nuclei decouple. 



prepared in a state with D — (Fig. [TJ:). In such a 
case, polarization stops and the dephasing time of the 
singlet-triplet qubit can be greatly extended. However, 
we have not found physical parameter regimes in which 
such states can be stably prepared. 

Model - The hyperfine coupling between a localized 
electron in dot d = l,r (for the left, right dot) and a nu- 
clear spin I kd at r kd , is given by g kd = a hi v \ip(r kd )\ 2 , 
where ip is the electron wavefunction, vq is the volume 
per nucleus, and ahf is a coupling constant. The homo- 
geneous limit is defined by g kd = g d for all fc. S and D 
are defined through the collective nuclear spin operators 
denoting the Overhauser fields in the left (L) and right 
(R) dots L = J2 k 9kilki and R = J2k 9k r Ikr such that 
S = (L + R)/2, D = (L - R)/2. 

For a double quantum dot with two electrons, we can 
write the Hamiltonian for the lowest energy (1,1) and 
(0, 2) electron states, where (n, m) indicates n (m) elec- 
trons in the left (right) dot. In this subspace the effective 
Hamiltonian for the electron and nuclear spins takes the 
form H = H e i + + H n , where 



E eX = 7e-B cx t • {at + s r ) + J(e)s e ■ s r 
i?hf = S ■ (se + s r ) + cos 0(e) D ■ (se — s r ) 



(1) 



k.d 



7n(-B Gx t + h kd ) ■ I, 



kd 



here sgM is the electron spin in the left(right) dot, 
7e {in) is the electron (nuclear) gyromagnetic ratio, 
where we consider spin 3/2 nuclei of a single species, 
-B C xt = B CKt z is the external magnetic field, cos 0(e) is 
the overlap of the adiabatic singlet state \s) with the 
(1,1) singlet state as a function of the detuning e be- 
tween the (1,1) and (0,2) singlet states, and J(e) is 
the splitting between \s) and |Tb) [2^|. The rms val- 
ues of the components of L, R in the infinite temper- 
ature ensemble arc Q d = (J2 k g kd I(I + l)/3) 1 / 2 . We 



define SI = ^/(ftf + tt 2 r )/2 w (10 ns)" 1 for typical few- 
electron double dot experiments, and work in units where 
Q = — 7 e = h = 1. In addition to the nuclear Zeeman 
energy we include a "noise" term h kd , representing the 
fluctuating, local magnetic field felt by a nuclear spin at 
site r kd , which could arise from e.g. nuclear dipole-dipole 
and electric quadrupole interactions. We estimate the 
scale of the fluctuations to be such that a typical nuclear 
spin dephases at a rate of 1-50 kHz (28l |. 

We find the nuclear spin evolution semiclassically by 
treating the nuclei and electrons as mean fields when solv- 
ing for the electron and nuclear dynamics, respectively. 
This scmiclassical approximation has been well studied 
in the context of central spin models and is generally 
reliable for extracting average quantities of high temper- 
ature, low polarization nuclear ensembles in dots with a 
large number of nuclei N (typically 10 6 [H]) 03, El- 

Neglecting H n , the nuclear spins evolve according to 
I kd = i[H h{ ,I kd ], giving equations of motion 



Ikd 



9kd 
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( (se + s r ) ± cos 6 [si — s T 



ikd) 



(2) 



where the top sign applies for d = I . We now replace 
(/fed) with I kd since we are treating the nuclear spins 
semiclassically. Consider a pulse cycle e{t) of duration 
T <C I /gkd ~ V^V /fl d . In a single cycle we can average 
over the fast evolution of the electrons to arrive at the 
coarse-grained equations (l2| 

i (+\ ~ Ikd{t + T) - hd{t) o f-t\ w t u\ r?\ 

Ikd{t) w — = gkd Pd(t) x Ikd(t), (3) 



T 



T+t dt' 



Pd(t)= I 



S r ) ± COS0 (S( — 8 r }], 



(4) 



where P d is a slowly-varying, effective Knight magnetic 
field felt by the nuclear spins. 

We now consider the class of pulse sequences employed 
in Refs. [T9I. [2lj] . in which the electronic system is initial- 
ized in \s) at large e and e is swept slowly through the 
|s)-|T_|_) resonance followed by a fast return to (0, 2) and 
reset of the electronic state via coupling to the leads. 
(Fig. QJi). This results in a build up of negative polar- 
ization. For simplicity, we work in the limit where the 
electron spin flip probability per cycle is small and cal- 
culate P d to lowest order in CI/ J, Cl/B eXt , CIT, and CI//3, 
where 1 = \ \dJ/dt\ \t=t r is the sweep rate at the reso- 
nance time t r , i.e., J[e(t r )] = B ext . 

To calculate (s d ) we work in the Heisenberg picture. 
Defining er™ = \T m ) (s\, we can write {sf — sf)/2 = {o-± — 
o-- 1 )/^ and (sf - s*)/2 = -(<r° + a n _)/2. Since £ cxt , 
J, /3 > CI, we can set (|T„) (T m \) = in (da™' /dt) to 
obtain the first order corrections to the electronic state: 

(a + ) = -iV2v(t)D z + iJ(t)(a° + ), (5) 
(a- 1 ) = -iv(t)D_ + i(J(t) + B cxt ) (a- 1 ) , (6) 
= iv(t)D+ + i(J(t) - B ext ) (a\_) , (7) 
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where v(t) = cos0(i)/\/2. Since J, B cxt > vQ, Eqs. [5] 
and [6] can be adiabatically eliminated. To find we 
formally integrate Eq. and perform a saddle point ex- 
pansion about the resonance time, assuming v(t) is con- 
stant in this region, to reduce it to a Landau-Zener prob- 
lem [2^] . From this solution we calculate the average ini- 
tial spin flip probability per cycle, Pfo = 2irv 2 (t r )il 2 / (3 2 . 
Putting these results into Eq. [4] gives 

P d = ±(r zxD ± -A DJ-A_D ± ), (8) 

where r = pfg/fl 2 T arises from the polarization process 
via T + , A = (2v 2 /J) c and A_ = (v 2 /(J + B ext )) c arise 
from electron-nuclear exchange processes via the To and 
T_ states, respectively, (-) c indicates an average taken 
over one cycle, and D± = D x x + D y y. 

Qualitatively, the effect of the To term is to polarize 
the nuclear spins, but it also saturates the polarization 
by driving the nuclear spins into a dark state, D±_ = 0. 
The Ao term drives the nuclear spins out of dark states, 
unless D z = as well. Without noise, states with D = 
are stationary during this DNP process; we refer to these 
as "zero states." 

Solving Eqs. [3J with P d given by Eq. [5J for an arbi- 
trary electron wave function is a challenging many-body 
problem. To help treat this problem, we have developed 
a new numerical method, which is formally equivalent 
to approximating the wave function by a unique set of 
M <C N coupling constants gkd, that well approximates 
the time evolution of L and R for a time that scales as 
M. A full description of this method, which was used 
in Fig. [2l along with a discussion of several higher order 
effects from finite magnetic field and adiabaticity, will be 
given elsewhere [j^ . 

Unequal dots - Our results that zero states are un- 
stable to the growth of large difference fields, in the 
presence of asymmetry in the size of the dots and nu- 
clear noise (H n ), can be shown analytically in the case 
of a simplified model. We assume homogeneous cou- 
pling and work in the high field, large J, limit where 
we can set Ao = A_ = in Pj. To treat the noise we 
first go into a frame rotating with the nuclear Larmor 
frequency, and assume can be rotated away. We 
further assume that the nuclear noise can be approxi- 
mated by a Gaussian, uncorrelated white noise spectrum, 
T^LW'w^'))" = 2nM - t')8 k k'5dd', where (•)„ are 
averages over the noise [30j. These local noise processes 
give rise to a mean decay of the collective nuclear spin 
variables L + (R+) and associated fluctuations J 7 ^), de- 
fined by {F d {t)F* d ,{t')) n = 2n 2 5 dd ,5{t - f ). As a result, 
Eqs. [3] and [51 including H n , give 

L+ = g e T L z (L + - R+)/2 - V L+ + ^Ti, (9) 

L z = -^T (Ll-R ± -L ± ), (10) 

and similarly for R. From Eq. [51 if we start in a zero 
state, Td will produce a fluctuation in D±, and the con- 
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FIG. 2: a) Long time limit of D z /S z as the relative hypcrfirie 
coupling in the two dots, R = g r /ge, is varied. The solid line is 
Eq.[T2]and the dashed line is ( 1 — i?) / ( 1 + 7?) , obtained from 
a heuristic model (see text). Circles are numerical results 
with statistical error bars after averaging over an ensemble 
of 1000 initial conditions, run out to t = 10 5 /g£ro ~ 1 s, 
using an approximation to a 2D Gaussian electron wavefunc- 
tion in terms of 100 coupling constants g^d with noise strength 
v/geT = 10~ 3 . b) Phase dia gram for identical dots for ei- 
ther saturation in dark states or the self-consistent growth of 
difference fields as the DNP pumping rate (vertical axis) and 
the Knight shift from |T_) (horizontal axis) are varied relative 
to the Knight shift from |Xb)- The dark grey shaded region 
is a numerical "crossover" regime where both effects occur 
depending on initial conditions and the dotted line is an an- 
alytic result from the simplified model of Eq. 1131 For typical 
polarization cycles A_/Ao ~ 1/4, but To/Ao « pfoB ext /Q 2 T 
can be tuned over a broad range. 

tribution to L z of the form —giToL^ results, in the long 
time limit, in L z <C —1 and similarly for R z . Thus, 
\L Z /L Z \ <C 1 and we can treat L Zi R z as static to find 
, and (L± ■ R±) n , which allow us to find 

the slow evolution of L z , R z . To lowest order in 1/S Z 
and 1 — R, where R = g r /gi, 

(D z ) n = -77 [(D z ) n - (1 - R) (S z ) n ] I (S z ) 2 n , (11) 

and (S z ) n = —y/rft. This growth of S z as t 1 / 2 is 
a result of our assumption of delta correlated nuclear 
noise. If we assume a finite correlation time r c such 
that {J 7 d{t)F* d {t')) n = Q2exp(- \t - t'\ /t c ) /t c , then for 
gT \S z \ < 1/t C) \S Z \ - i 1 / 2 , but eventually \S Z \ - t 1 / 3 . 
Integrating Eq.fTTI gives D z /S z — > (l—R)/2. For general 
R we find, in the long time limit, 



S s 2R+^4R 2 + (l-R) 4 ' 

Fig. [5^ shows good agreement between these results and 
numerics for an inhomogeneous Gaussian wavefunction. 

Identical dots - For identical dots the previous argu- 
ments are no longer valid. Fig. [2b, however, shows the 
results of numerical simulations [29j that demonstrate the 
existence of a parameter regime for which there is self- 
consistent growth of D z even for identical dots. Simula- 
tions were performed at each set of parameters by taking 



4 



20 different initially polarized nuclear spin configurations 
with S g = -10, D z = -2, ri/giAo between 10~ 2 - 10~ 4 , 
and a 2D Gaussian electron wavefunction approximated 
with 400 values of gkd- We determined which parame- 
ter values had (D z ) growing after t = 10 3 /.g;Ao. For 
To/Ao > 1/2, no self-consistent growth of D z appears, 
and the system approaches a dark state. For smaller 
r /A and for moderate A_/A , continued growth of 
D z is observed. We find a similar boundary for unequal 
dots provided |1 - R\ < 0.05. 

This phase diagram for identical dots can be verified 
analytically in a simplified model, where the hypcrfine 
coupling in each dot takes two values {g\ 3> g 2 ) on 
two groups of spins of similar size. We assume initially 
—g 2 S z ^> g\ \D Z \ ^> 1 ^> D±_ with the polarization 
mostly in the strongly coupled spins. To lowest order 
in g 2 /gi, r//g 2 D z , giD z /g 2 S z , and Dj_/D z , we find [2^ 

(£>,)» ex (Tg + A 2 „ - A A_)( 5l (D z ) n /g 2 (S z ) 7 f. (13) 

Growth of D z requires nonzero D±, but, as we show 
below, for large polarization and weak noise D± ~ 
D z /S z , which implies that the growth D z must occur 
self-consistently to prevent saturation. This is illustrated 
by Eq. 1131 where the continued growth of D z is entirely 
determined by the sign of Tg + A 2 . — AoA_. For large 
To or strong DNP pumping, the sign is positive, satu- 
ration effects dominate, large difference fields are unsta- 
ble and the system eventually reaches a dark state. For 
smaller Tq, the sign is negative and coherent evolution 
arising from interactions with iTb,-) allows D z to con- 
tinue growing and D± remains finite. Fig. [2Jd shows rea- 
sonable agreement between this predicted boundary and 
our numerical results. 

We now address the stability of the zero states in the 
absence of nuclear noise. For identical dots, in the ho- 
mogeneous limit, Eqs. [3] and [8] give 

D+ = gi{&- - iT )S z D + ~giA Q D z S+, (14) 
D z = g [(A_ - ir )D + S_ - c.c] /2i. (15) 

Near a zero state S is constant since S ~ 0(D 2 ). The po- 
larization, gToS z , acts as a damping term for £)+; conse- 
quently, for S z « -1, D+ -> [A 5+/(A_ - iT )]D z /S z . 
Together with Eq. [TS] this implies D z = 0. Thus the 
stability matrix, <9-D m /<9_D„|.d = o, has two negative eigen- 
values and one zero eigenvalue. Due to this zero eigen- 
value, we expect the stability of a zero state to be highly 
sensitive to external perturbations. We find that inhomo- 
geneous hyperfine coupling, multiple nuclear species, the 
hybridization of \s) and |T ) as discussed in Refs. [25|.|26|. 
and additional higher order corrections to P c i in 1/B cxt 
do not, however, break this zero eigenvalue. In the ab- 
sence of noise, we find numerically that for some param- 
eters a large fraction of initial conditions result in the 
system spending a long time near a zero state; however, 



when we include nuclear noise or higher order corrections 
in the inverse sweep rate, for example, zero states become 
repulsive on a long time scale [2^ . Throughout this work 
we have mostly neglected nuclear spin diffusion [30I ] and 
spin-orbit coupling [18[ , both of which could potentially 
affect DNP and, in particular, the stability of zero states. 
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[1 
[2 
[3] 
[4] 

[5; 

[6 
[7; 

p: 

[9 
[10 

[11 
[12 

[13 

[14; 

[15 

[16 

[17; 

[18 
[19 
[20 
[21 
[22 
[23 

[24 

[25 
[26 

[27 

[28 
[29 
[30 



G. Yusa et al, Nature 434, 1001 (2005). 

D. C. Dixon et al., Phys. Rev. B 56, 4743 (1997). 

G. Salis et al., Phys. Rev. Lett. 86, 2677 (2001). 

K. Ono and S. Tarucha, Phys. Rev. Lett. 92, 256803 

(2004). 

F. H. L. Koppens, K. C. Nowack, and L. M. K. Vander- 
sypen, Phys. Rev. Lett. 100, 236802 (2008). 

A. S. Bracker et al., Phys. Rev. Lett. 94, 047402 (2005). 

C. W. Lai et al., Phys. Rev. Lett. 96, 167403 (2006). 
W. A. Coish and J. Baugh, Phys. Status Solidi B 246, 
2203 (2009). 

R. Hanson et al., Rev. Mod. Phys. 79, 1217 (2007). 

G. Chen, D. L. Bergman, and L. Balents, Phys. Rev. B 
76, 045312 (2007). 

K. A. Al-Hassanieh et al., Phys. Rev. Lett. 97, 037204 

(2006) . 

H. Christ, J. I. Cirac, and G. Giedke, Phys. Rev. B 75, 
155324 (2007). 

W. M. Witzel and S. D. Sarma, Phys. Rev. B 77, 165319 

(2008) . 

W. Yao, R.-B. Liu, and L. J. Sham, Phys. Rev. B 74, 
195301 (2006). 

W. A. Coish and D. Loss, Phys. Rev. B 70, 195340 
(2004). 

D. Klauser, W. A. Coish, and D. Loss, Phys. Rev. B 78, 
205301 (2008). 

J. R. Petta, H. Lu, and A. C. Gossard, Science 327, 669 
(2010). 

M. Rudner et al. (2009), arXiv:0909.0060. 

S. Foletti et al., Nat. Phys. 5, 903 (2009). 

J. R. Petta et al., Phys. Rev. Lett. 100, 067601 (2008). 

D. J. Reilly et al., Science 321, 817 (2008). 

S. Foletti et al. (2008), arXiv:0801.3613. 

G. Ramon and X. Hu, Phys. Rev. B 75, 16130RR) 

(2007) . 

H. Ribeiro and G. Burkard, Phys. Rev. Lett. 102, 216802 

(2009) . 

W. Yao (2009), arXiv:0905.2460. 

M. Stopa, J. J. Krich, and A. Yacoby, Phys. Rev. B 81, 
041304(R) (2010). 

J. M. Taylor, A. Imamoglu, and M. D. Lukin, Phys. Rev. 
Lett. 91, 246802 (2003). 

J. M. Taylor et al., Phys. Rev. B 76, 035315 (2007). 
M. Gullans et al., in preparation. 

D. J. Reilly et al., Phys. Rev. Lett. 101, 236803 (2008). 



